Octet baryon magnetic moments from QCD sum rules 
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A comprehensive study is made for the magnetic moments of octet baryons in the method of 
QCD sum rules. A complete set of QCD sum rules is derived using the external field method 
and generalized interpolating fields. For each member, three sum rules are constructed from three 
independent tensor structures. They are analyzed in conjunction with the corresponding mass sum 
rules. The performance of each of the sum rules is examined using the criteria of OPE convergence 
and ground-state dominance, along with the role of the transitions in intermediate states. Individual 
contributions from the u, d and s quarks are isolated and their implications in the underlying 
dynamics are explored. Valid sum rules are identified and their predictions are obtained. The 
results are compared with experiment and previous calculations. 

PACS numbers: 13.40.Em, 12.38.-t, 12.38.Lg, ll.55.Hx, 14.20.Gk, 14.20.Jn 



I. INTRODUCTION 

The QCD sum rule method is a nonperturbative ana- 
lytic formalism firmly entrenched in QCD with minimal 
modeling. The field remains active judging by the 3000 
and growing references to the seminal paper of Shifman, 
Vainshtein and Zakharov [l[ that introduced the method. 
The approach provides a general way of linking hadron 
phenomenology with the interactions of quarks and glu- 
ons via only a few parameters: the QCD vacuum con- 
densates and susceptibilities. The studies give an unique 
perspective on how the properties of hadrons arise from 
nonperturbative interactions in the QCD vacuum and 
how QCD works in this context. It has been success- 
fully applied in almost every aspect of strong-interaction 
physics. 

Calculations of the magnetic moment were carried out 
soon after the method was introduced for the proton, 
neutron 0, H| and hyperons [1] in the external field 
method. In this method, a static magnetic field is in- 
troduced that couples to the quarks and polarizes the 
QCD vacuum. Magnetic moments can be extracted from 
the linear response to this field. The results of the studies 
validated the external field method as a way of probing 
hadron properties other than the mass, such as magnetic 
moments, form factors, axial charge, isospin- breakings. 
Later, a more systematic study was made for the mag- 
netic moments of octet baryons 0, H, 0, @| • Calculations 
were also carried out for decuplet baryons [H, [l(| HH, E3 
and the rho meson [TI|. There are other studies of 
magnetic moments using the light-cone QCD sum rule 
method [3, [TH, [l7| which will not discuss here. 

In this work, we carry out a comprehensive, indepen- 
dent calculation of the magnetic moments of the octet 
baryons in the external field method. It can be consid- 
ered as an update over the previous calculations [U, H 0] 
which were done more than 20 years ago. There are a 
number of things we do differently. First, we employ 
generalized interpolating fields which allow us to use the 
optimal mixing of interpolating fields to achieve the best 
match. Second, we derive a new, complete set of QCD 



sum rules at all three tensor structures and analyze all 
of them. The previous sum rules, which were mostly 
limited to one of the tensor structures, correspond to 
a special case of the mixing in our sum rules. In this 
way, we provide an independent check of the previous 
sum rules. Third, we perform a Monte-Carlo analysis 
which has become standard nowadays. The advantage 
of such an analysis is explained later. Fourth, we use a 
different procedure to extract the magnetic moments and 
to treat the transition terms in the intermediate states. 
Our results show that these transitions cannot be simply 
ignored. Fifth, we isolate the individual quark contribu- 
tions to the magnetic moments and discuss their impli- 
cations in the underlying quark-gluon dynamics in the 
baryons. 

The paper is organized as follows. In Section |TT1 the 
method of QCD sum rules is introduced. Using the in- 
terpolating fields for the octet baryons, master formula 
are calculated. Then both the phenomenological repre- 
sentation and QCD side are derived. Section ILTIl will list 
the sum rules we derived for the octet baryon family, fol- 
lowed by the analysis to extract the magnetic moments in 
Section HVl Section IVl summarizes the results and gives 
a comparison of our results with experiment and previ- 
ous calculation, followed by an in-depth discussion of our 
findings. Our conclusions are given in Section IVH 



II. METHOD 

The starting point is the time-ordered correlation func- 
tion in the QCD vacuum in the presence of a constant 
background electromagnetic field F M „: 



n(p) = i d 4 x e ip x (0 | T{ r)(x) fj(0) } | 0)j 



(1) 



The QCD sum rule approach is to evaluate this corre- 
lation at two different levels. On the quark level, it 
describes a hadron as quarks and gluons interacting in 
the QCD vacuum. On the phenomenological level, It 
is saturated by a tower of hadronic intermediate states 
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with the same quantum numbers. This way, a connec- 
tion can be established between a description in terms of 
hadronic degrees of freedom and one based on the under- 
lying quark and gluon degrees of freedom governed by 
QCD. Here ij is the interpolating field (or hadron cur- 
rent) with the quantum numbers of the hadron under 
consideration. The subscript F means that the correla- 
tion function is to be evaluated with an electromagnetic 
interaction term added to the QCD Lagrangian: 

d = -A^J", (2) 

where A^ is the external electromagnetic potential and 
= BqQ'yPq is the quark electromagnetic current. 
Since the external field can be made arbitrarily small, 
one can expand the correlation function 

n(p)^u^( P ) + n^(p) + --- , (3) 

where IT^ ) (p) is the correlation function in the absence 
of the field, and gives rise to the mass sum rules of the 
baryons. The magnetic moments will be extracted from 
the QCD sum rules obtained from the linear response 
function H^(p). 

The action of the external electromagnetic field is two- 
fold: it couples directly to the quarks in the baryon in- 
terpolating fields, and it also polarizes the QCD vacuum. 
The latter can be described by introducing new parame- 
ters called vacuum susceptibilities. 

The interpolating field is constructed from quark fields 
with the quantum number of baryon under consideration 



and it is not unique. We consider a linear combination 
of the two standard local interpolating fields. They read 
for the baryon octet family: 



r]P(uud) = ~2e abc [(u aT Cj 5 d b )u c + f3(u aT Cd b )j 5 u% 
rj n (ddu) = -2e abc [(d aT Cj 5 u b )d c + (3{d aT Cu b ) lb d% 

n A (uds) = -2^e abc [2(u aT Cj 5 d b )s c + {u aT C l5 s b )d c 

-(d aT Cj 5 s b )u c + {3{2{u aT Cd b ) l5 s c + {u aT Cs b ) l5 d c 
-(d aT Cs b ) l5 u c )}, 

rf {dds) = -2e abc [{d aT C~f 5 s b )d c + f3{d aT Cs b )-f 5 d% 
v *°(uds) = -V2e abc [(u aT Cj 5 s b )d c + {d aT C lb s b )u c 
+p{(u aT Cs b ) l5 d c + {d aT Cs b )j 5 u c )l 
T]^ + (uus) = -2e abc [(u aT C lb s b )u c + (3(u aT Cs b )j 5 u c ], 
rf~ {ssd) = -2e abc [{s aT C l5 d b )s c + (3{s aT Cd b ) l5 s% 
rf\ssu) = -2t ab %s aT C^u b )s c + /3(s aT Cu b )~f 5 s c }. 

(4) 

Here u and d are up-quark and down-quark field oper- 
ators, C is the charge conjugation operator, the super- 
script T means transpose, and e a & c makes it color-singlet. 
The normalization factors are chosen so that correlation 
functions of these interpolating fields coincide with each 
other under SU(3)-flavor symmetry. The real parameter 
[3 allows for the mixture of the two independent currents. 
The choice advocated by Ioffe [l8| and often used in QCD 
sum rules studies corresponds to (3 — — 1. We will take 
advantage of this freedom to achieve optimal matching 
in the sum rule analysis. 



I 

A. Phenomenological Representation 

We start with the structure of the two-point correlation function in the presence of the electromagnetic vertex to 
first order 

U(p) = i J d 4 xe^ x < Q\ V (x)[-i J d 4 yA^y)J^(y)}fj(0)\0 > ■ (5) 

Inserting two complete sets of physical intermediate states, we restrict our attention only to the positive energy ones 
and write 

U{p) = J d 4 xd 4 y-^-^ Y^J2 k ,2_ M l^ lE k 2_ M l_ is 
e l P' x A^{y) < 0\r](x)\N'k's' >< N'k' s'\J»(y)\Nks >< 7Vfcs|77(0)|0 > . 
We can use the translation invariance to express rj{x) in terms of i](0) 

< Olrjix^N'k's 1 >=< 0\r)(0)\N'k's' > e~ ik ' x . (7) 

The interpolating field excites (or annihilates) the ground state as well as the excited states of the baryon from the 
QCD vacuum. The ability to do so is described by a phenomenological parameter Ajv (called current coupling or pole 
residue), defined by the overlap for the ground state 

< 0|?7(0)|7Vfcs >= X N u(k, s), (8) 

where u is the Dirac spinor. 

Translation invariance on J^(y) gives 



< N'k's'\ J»{y)\Nks >= e l ™ < N'k's'\ J>*(0)\Nk8 >, 



(9) 
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where q — k! — k is the momentum transfer and Q 2 = —q 2 . 

The matrix element of the electromagnetic current has the general form 

< k's'\J^0)\ks >= u(k', s'^F^Q 2 )^ + F 2 (Q 2 )ia^^}u(k, s), (10) 

where the Dirac form factors F\ and F 2 are related to the Sachs form factors by 



G £ (Q 2 ) = F 1 (Q 2 )- I5 ^F 2 (Q 2 ) 
G M (Q 2 )^F 1 (Q 2 ) + F 2 (Q 2 ). 



(2M») ar2 W ) III) 



At Q 2 = 0, Fi(0) = 1, F 2 (0) = [i a which is the anomalous magnetic moment, and Gm(0) = i*i(0) + ^2(0) = /1 which 
is the total magnetic moment. 

Writing out explicitly only the contribution of the ground-state nucleon and denoting the excited state contribution 
by ESC, we have 

n(p) - -\%Jd l xd i y ]$ H $ i [k* -M%- ie]^[k 2 -M 2 N - ie]" 1 ^)^-*')* 

e*™ £ *W> s')[^i(0 2 )7 M + F 2 {Q 2 )io^^\ £ u(fc, «)u(fc, a) + ESC. ^ 

s' s 

The spin sums are of the form 

£w(fc,s)u(fc,s) = k + M N . ( 13 ) 

s 

QCD sum rule calculations are most conveniently done in the fixed-point gauge. For electromagnetic field, it is defined 
by x^A^x) = 0. In this gauge, the electromagnetic potential is given by 

Mv) = -\f^V V - (14) 
Changing variables from k to q = k' — k, then d 4 k = —d 4 q, we have 

n(p) = =$LF„, fd^y^^W* - M 2 N - ze]^[{q - k') 2 - M% - l e]^ 
xe^-^i-i^e^ip + M 7V )[F 1 (Q 2 ) 7ai + F 2 {Q 2 )iu» v ^]{k' - q + M N ) + ESC. 

Integrating over x, we get a delta function 

/ d4 *ji^e i{p - k ' )x =5 4 (P-k'). (16) 

Integrating d jdq v by parts, then doing J d 4 y, we can get another delta function 6(q). Since we have a 6(q), when doing 
d/dq v only terms linear in q v contribute. We have i*i(<3 2 )| 9 =o = 1, F 2 (Q 2 )\ q =o = fi a = /i— 1, and d / dq„Fi t2 (Q 2 )\ q=0 = 
0, so derivatives of the structure functions do not enter. Finally, we arrive at 

n(p) = ^\ 2 N F^[p 2 -M 2 N -ie]- 2 {p + M N ){^^a^(p + M N )+i<j^ . 
-(p"7 v -p v ^){p + M N )[p 2 - M 2 - is}- 1 } + ESC. { ' 

Examination of its tensor structure reveals that its has 3 independent combinations: F ,iV (pa llv + cr^p), F^ ,v i(p tl j v — 
Pvlp)p and F^Opv. The momentum-space correlation function in the above equation can be written in terms of 
these three structures 



H(p) = - ^-fc P {K1(2M^) + ^{p 2 - M 2 N ) + m \pa^ + a^p\ 

+^ j ^[i(P tt r-p"in]p}+Esc, 

where we have used the following identities 



(p + M N )a»" = M N a» v + \{pa^ v + a^p) + i(p»-y" - p u ^), 

(p + M N )cr^(p + M N ) = a^(p 2 + M 2 N ) + M N {pa^ + u^p) + 2i{p^ v - p v ^)p, (19) 
(p + M w )(p" 7 " - p" 7 ")(p + M N ) = -(p 2 - M 2 N )(p»Y - p v ^). 
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FIG. 1: The three kinds of contributions to the spectral function in the presence of an external field: ground state, transitions 
between ground state and excited states, and pure excited states 



Next step is to perform the Borel transform defined by 



B[f(p 2 )} = 



lim 

-p 2 ,n — >oo 
-p 2 /n=M' 



(20) 



Upon Borel transform the ground state takes the form 



B[IL(p)] 



iM 2 



(21) 



where M is the Borel mass, not to be confused with the nucleon mass Mm- 

Here we must treat the excited states with care. For a generic invariant function, the pole structure can be written 

as 



r- 



(p 2 -M 2 N ) 2 



E 

N* 



(p 2 



r 2 

N ^ N* 

-M 2 N )(p 2 -M 2 ,) 



+ E 



N* 



(pi-M 2 ,) 2 



(22) 



where Cn^n, C'n<^n* and Cjv»<-»at» are constants. The first term is the ground state pole which contains the desired 
magnetic moment fj,N- The second term represents the non-diagonal transitions between the ground state and the 
excited states caused by the external field. The third term is pure excited state contributions. These different 
contributions can be represented by the diagrams in Fig. [TJ Upon Borel transform, it takes the form 



-Ml/M 2 , -M 2 N /M 2 



M 2 



E 



Mi 



'-M% ( 



_ e -{M 2 N ,-M 2 N )/M-- 



E 



M 2 



-M 2 N ,/M 2 



(23) 



The important point is that the transitions give rise to a contribution that is not exponentially suppressed relative 
to the ground state. This is a general feature of the external-field technique. The strength of such transitions at 
each structure is a priori unknown and is an additional source of contamination in the determination of the magnetic 
moment /ijv- The standard treatment of the transitions is to approximate the quantity in the square brackets by a 
constant, which is to be extracted from the sum rule along with the ground state property of interest. Inclusion of such 
contributions is necessary for the correct extraction of the magnetic moments. The pure excited state contributions 
are exponentially suppressed relative to the ground state and can be modeled in the usual way by introducing a 
continuum model and threshold parameter. 



B. Calculation of the QCD Side 

We start by contracting out the quark pairs in Eq. {T]) using Wick's theorem, resulting in the so-called master 
formula in terms of quark propagators. The master formula for the proton (with uud quark content) is 

(n | T{r) N (x) fj N (0) | 0) = _4e Qbc e Q ' bV { 
+Sf l5 CSf T C l5 S b u b ' + SfTr(CSf T C l5 S b u b ' l5 ) 
+Pl 5 Sf l5 CSf T CS b u b ' + (3 l5 SfTr(CS c u c ' T CS bb ' l5 ) 
+(3Sfcsf T C l5 S b u b ' l5 + [iSf lb TT{CS c u c ' T C l5 S bb ') 

+(3 2 l5 SfCSf T CSf l5 + p 2 l5 Sf l5 TT(CSf T CSt b ') }. (24) 
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The master formula for neutron (with ddu quark content) can be obtained by exchanging the d quark with a u quark 
from Eq. (??). By replacing the d quark with a s quark, one can get the master formula for S + (with uus quark 
content). While the master formula for £~ (with uus quark content) can be obtained by replacing the u quark with 
a d quark from E + master formula. By exchanging the u quarks and s quarks in E + master formula, the master 
formula for S° (with ssu quark content) can be obtained. Likewise, by replacing the u quark with a d quark, one can 
get the master formula for S + (with ssd quark content). 

The master formulae for E° and A (uds quark content) have a more complicated structure. They can be written 
in a combined way as 

(n\T{ T 1 (x)fj(0)}\n) = -fe abc e a ' b ' c '{ 

fif 2 sfj 5 cs c u c ' T c l5 s bb ' - hf 3 sf l5 csT' T c l5 s bb ' + hhsf^{cs c u c ' T c lb sf l5 ) 

+.f 2 fiSf l5 CS c u c ' T C l5 Sf + hhSf^CSf T C lb sf + f 2 f 2 S a /Tr(CSf T C l5 Si b ' l5 ) 

-f 3 hsf l5 csf T c l5 sf + f 3 f 2 sf l5 csf T c l5 s b / + hhsf^{csf T c lb sf lb ) 

+hh(3SfCSf T C l5 S b / l5 - f 1 f 6 (3SfCSf T C l5 St b ' l5 + hf 4 psf l5 Ti{CS c d c ' T C l5 S bb ') 
+hf4[3SfCS c u c ' T CToS bb ' l5 + f 2 f e (3SfCSf T C To S b u b ' l5 + f 2 f 5 (3Sf l5 Tr(CSf T C l5 S b u b ') 
-f 3 f4(3Sfcsr' T C l5 S bb ' l5 + f 3 f 5 (3Sfcsf T C l5 S bb ' l5 + hhPSf l5 Tr(CSf T C l5 S bb ') 
+hhPl 5 S a s a 'lsCS c u c ' T CS bb ' - f 4 f 3 (3 l5 sf l5 CS c / T CSf + A/i/375Sr'Tr(C^ c ' T CSf 75) 
+hh(3l 5 sf l5 CSf T CS bb ' + f 5 f 3 (3 l5 sf l5 CSf T CSt b ' + hf 2 p lb S a /^{CSf T CS bb ' lb ) 
-hfi(3l 5 Sf l5 CSf T CS bb ' + hf 2 (3 l5 sf l5 CSf T CS bb ' + hf 3 p l5 SfTr(CSf T CS bb ' l5 ) 
+f4h(3 2 l 5 SfCS c u c ' T CS bb ' l5 - f 4 f 6 (3 2 l5 sfcS c / T CSi b ' l5 + f 4 f 4 (3 2 l5 Sf l5 Tr(CSf T CS bb ') 

+hf 4 (3 2 l5 sfcs c u c ' T cs bb ' l5 + hM 2 l5 sfcsf T csi b ' l5 + hM 2 l5 sf l5 Ti{CSf T CS bb ') 

-hh^S a u a 'CST' T CS bb ' lb + f 6 f 5 p 2 -f 5 sfcsf T CS bb ' l5 + Af 6 f 6 ^ 5 sf l5 Tv(CSf T CS bb ')}, (25) 

where the various factors are: for £°, / = 2/3, fi — 0, f 2 — 1, f 3 = 1, f 4 = 0, /s = 1, /g = 1; and for A, / = 2, /1 = 
2,/ 2 = l I /3 = -l,/4 = 2,/ 5 = l,/ 6 = -l. 
In the above equations, 

S a q b (x,0;F) = (0\T{q a (x)cf'(0)}\0) F , q = u,d,s, (26) 

is the fully interacting quark propagator in the presence of the electromagnetic field. To first order in and m q 
(assume m u — ma = 0, m s ^ 0), and order x 4 , it is given by the operator product expansion (OPE) 0, H, 

+ ^9c.o ■ Gq)x 2 6 ab - ^(? M • Gq)xx 2 5 ab - ^{qq)( g 2 c G 2 )x 4 5 ab 
1 xa^ + a^x fX n \ ab 1 * 2ri2s xa a P +a a t 3 x f \ n 



-> ( T J +48 32^ G 

lb -, / \n \ Q>b 

3 2 2 10 



+ f> c * • Gq) + ^x) (£) ^ + ^^±^S« b 

-^L x (qq)F aP o^5 ab + l ^l x (qq)F a0 (i^ + a^x) 5 ab 
+^q)F a0 (x 2 *^ - 2x p W) 5 ab 
+^(qq)F a0 [x 2 ( K + 0<J afj ~ x p x' } (2k - 6 ab 

\n \ ab 



~(qq) ( KF a p - -& a p^F^ ) ( — ) + higher order terms. (27) 
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We use the convention g 0123 = +1 in this work. 



In addition to the standard vacuum condensates, the 
vacuum susceptibilities induced by the external field are 
defined by 

{qa» v q)F = e q x{qq) F ^v, 

{qg c G flu q) F = e q n(qq)F^, (28) 
{qg^tivpxG^^^F = ieqZiqqjF^. 

Note that x has the dimension of GeV -2 , while k and £ 
are dimensionless. 

With the above elements in hand, it is straightforward 
to evaluate the correlation function by substituting the 
quark propagator into the various master formulae. We 
keep terms to first order in the external field and in the 
strange quark mass. Terms up to dimension 8 are consid- 
ered. The algebra is extremely tedious. Each term in the 
master formula is a product of three copies of the quark 
propagator. There are hundreds of such terms over vari- 
ous color permutations. The calculation can be organized 
by diagrams (similar to Feynmann diagrams) in Fig. [2] 
and Fig. [3) Note that each diagram is only generic and all 
possible color permutations are understood. The QCD 
side has the same tensor structure as the phenomenolog- 
ical side and the results can be organized according to 
the same 3 independent structures. 



III. QCD SUM RULES 

Once we have both the QCD side (LHS) and the phe- 
nomenological side (RHS), we can derive the sum rules by 
matching both sides. Since there are three independent 
tensor structures, three sum rules can be constructed. 
We denote these tensor structures by the following short- 
hand notation 

WEi = F^ipa^ + a^p), (29) 
WOi = F^a^, (30) 
W0 2 = F^iip^-p^p (31) 

The sum rule from WEi involves only dimension-even 
condensates, so we call this sum rule chiral-even. The 
sum rule from both WOi and WO2 involves only 
dimension-odd condensates, so we call them chiral-odd. 
Note that in previous works [H, Q the dimension of the 
tensor structures, rather than the dimension of the con- 
densates, was used to refer to the sum rules. The two 
names are opposite. 

Now we are ready to collect all of the QCD sum rules. 
At the structure WEi, all of the sum rules can be ex- 
pressed in the following form 



Cl L- 4 / 9 E 2 (w)M 4 + c 2 m s xaL- w / 27 E 1 {w)M 2 + c 3X a 2 L-^ 27 E {w) + Ci bL- 4 / 9 E Q (w) + (c 5 + c 6 )m s aL~^ 9 E (w) 
+ (c 7 + c 8 )a 2 i 4 / 9 ^ + cgxm^L-W^ + Clomsm 2 ai -26/27_^ 

= -\ 2 N [^+A]e- M »' M \ 

, (32) 



where the coefficients differ from member to member. 
The quark condensate, gluon condensate, and the mixed 
condensate are 

a = — (2tt) 2 (uu), b=(g 2 G 2 ), (ug c a ■ Gu) = —to 2 , (uu) . 

(33) 

The quark charge factors e q are given in units of electric 
charge 



e u = 2/3, e d = -1/3, 



-1/3. 



(34) 



Note that we choose to keep the quark charge factors 
explicit in the sum rules. The advantage is that it can 
facilitate the study of individual quark contribution to 
the magnetic moments. The parameters / and (f> account 
for the flavor-symmetry breaking of the strange quark in 
the condensates and susceptibilities: 



/ = 



(ss) (sg c o ■ Gs) 
(uu) (ug c a ■ Gu) 



Xs K S £s rt r\ 

= — = — = 7"- (35) 



The anomalous dimension corrections of the interpolat- 
ing fields and the various operators are taken into account 



in the leading logarithmic approximation via the factor 



L 1 = 





7 


~\n{M 2 /k 2 QCD ) 


7 


a s (M 2 )_ 









(36) 



where [i = 500 MeV is the renormalization scale and 
Aqcd is the QCD scale parameter. As usual, the pure 
excited state contributions are modeled using terms on 
the OPE side surviving M 2 — > 00 under the assumption 
of duality, and are represented by the factors 



E n (w) = 1 



/m 2 v (w 2 /M 2 r 



(37) 



where w is an effective continuum threshold and it is in 
principle different for different sum rules and we will treat 
it as a free parameter in the analysis. 

Also, X% is the rescaled current coupling Xjy = A 2 /4. 

We only need to carry out four separate calculations 
for E+(uus), S°(uds), S+(ssd) and A(uds). The QCD 




FIG. 3: Diagrams considered for the strange quark mass corrections to the octet baryon magnetic moments. 
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sum rules for other members can be obtained from them by making appropriate substitutions specified below. 



Cl = i[(l + /3) 2 c s - 2(3 + 213 + 3/3 2 )c M ]; 

C2 = 4(1 + /?)[(-! + 3/3)c s / s </> - 2(1 + 3/3)c„]; 

c 3 = ti(-1 + /3)[(1 + /3)c s / 5 - 2c„(l -/? + /,+ /?/.)]; 

c 4 = -sIj [(3 - 2/3 + 3/3 2 )c s + 2(11 + 6/3 + ll/3 2 ) Ctt ]; 

c 5 = £[(3 - 18/3 + 3/3 2 - 2f s + 2/3 2 / s )c s + (-5 + 2/3 + 27/3 2 - 9/ s + 60 f a + 27/3 2 / s )c u ]; 

c 6 = + + P) c sf*<P( &K - + 2c «( 4 « + 6/3k + 0]; 

c 7 = ^(-1 + /3)[c s (-3 + 3/3 + 2/ s + 2/?/ s ) + 2c„(-l + /3 + 8/ s + 8/3/ s )]; 
c 8 = -^(-1 + /3)[(1 + /3)c s / s 0(8« - 

+c u (-2«: + 2/3 K + 10/ sK + 10(3f s K + £ - ft + ,U + Pf,Z)]; 
c 9 = -mi( 7 ^ + P)c s fs^ - c u (3 - 3/3 + f s + /3/ s )]; 
cio = -i[(l - 14/3 + /3 2 )c s + (-23 - 2/3 + 13/3 2 )c„(l + /,)]. 



ci = ^[(l + /3) 2 c„-2(3 + 2/3 + 3/3 2 )c s ] ; 

c 2 = 4(1 + /?)[(-! + 3/3)c u - 2(1 + 3/3)c s / s 0]; 

c 3 = i(-l + /?)[(1 + /3)c„/ s - 2c s / a 0(/ 8 - /3/ s + 1 + /3)]; 

c 4 = -skP - 2/3 + 3/3 2 )c u + 2(11 + 6/3 + ll/3 2 )c s ]; 

c 5 = 2i[((3 - 18/3 + 3/3 2 )/ s - 2 + 2/3 2 )c s + ((-5 + 2/3 + 27/3 2 )/ s - 9 + 6/3 + 27/3 2 )c s 

c 6 = -^(-! + /?)[(! + 0)c„(8k - + 2c s /^(4« + 6/3« + 0]; 

c 7 = ^(-1 + /3)[c„.f s ((-3 + 3/3)/ s + 2 + 2/3) + 2c s / s ((-l + /?)/, + 8 + 8/3)]; 

c 8 = -^(-1 + /3)[(1 + /3)c„/ s (8« - 

+c s f s cf>(-2 K f s + 2/3k/ s + 10k + 10/3« + £ - 0f s £ + f s £ + ft)}; 

c * = -2S8 1( 7 (! + - c -/^(( 3 - S/ 3 )-^ + 1 + 0)Y> 

cio = -i[(l - 14/3 + /3 2 )c, u + (-23 - 2(3 + 13/3 2 )c s (l + /,)]. 



c i = 'Mi' 1 + Pf( c d + c u) + (13 + 10/3 + 13/3 2 )c s ); 

c 2 = -Wl + /3)((-l + /3)(c d + c) - (1 + 50)c a f 8 4>); 

c 3 = + /?)((<* + c„)(l - /3 + / s + 5/3/ s ) - 2(5 + /3)c s / s </>); 

c 4 = - nk(( 17 + 2 ^ + 17 /? 2 )( c rf + c «) + (41 + 26/3 + 41/3 2 )c s ); 
c 5 = I l a (c s (-78 + /3(84- 16/ s )+/3 2 (210-4/ s ) + 20/ s ) 

+ (-1 + 0)(c d + c)(7 + 3/ s + /3(41 + 33/,))); 
c 6 = -^(-1 + /3)(c s / s 0(8(l + 2/3)k + (5 + /3)0 + (cj + c„)((20 + 22/3)k - (1 + 2/3)0); 
c 7 = ^(-1 + /3)(c s (39 + 33/3 - 10/ s - 20f s ) + (c d + c u )(l - (3 + 4/ s + 20/3/,)); 
c s = -252 (- 1 + /3)(2c s / s ^(8(l + 20)k + (5 + /?)£) 

+(c d + c„)(2(19 + 17/3 + / s + 5/3/ s )k - (1 - + f s + 5/3/ s )£)); 
c 9 = (-1 + /3)((c d + c M )(5 + 3/3 + 3/ s + 7/3/ s ) - 2(5 + 0)c s f s <j>); 
c w = ^((62 - 4/3 - 34/3 2 )c s - (-7 - 10/3 + 5/3 2 )(c d + c„)(l + /,)). 



At the structure WOi, the sum rules can be expressed in the following form 



Cl m s L- 8 / 9 E 2 (w)M 4 + c 2 xaL- 1 ^ 27 E 2 (w)M 4 + (c 3 +c 4 )a£a(w)M 2 + c 5 m sX a 2 L- 16 / 27 E (w) 



For £+ at WE i: 



For S° at WE i: 



For A at WEi: 




(41) 



For £+ at WOi: 

ci = -±(-l + /3)((l + /3)c 5 -2c u ); 

c 2 = -£(-1 + /?)[(- 1 + /3)c s / s + 18(1 + /3)c„]; 

c 3 = £(-1 + /?)[c s (-6 - 6/3 - / s + /3/ s ) - 6c u (2 + 2/3 - / s + /?/,)]; 

c 4 - i(-l + /3)[(-l + /?)c a / a 0(14/« - 130 - 18(1 + /3)c„(2k + £)]; 

c 5 = -^[(3 + 2/3 + 3/3 2 )c 5 / s </> - (3 + 4/3 + 9/3 2 )c M (l + /.)]; 

c 6 = i (-1 + /3)[(-l + /3)c s / s </> - 18(1 + /3)c„]; 

c 7 - -i[c 8 (-3 + 3/3 2 + 5/ s - 2/3/ s + 5/3 2 / s ) - c„(3 + 4/3 + 9/3 2 - 3/ s + 10/3/ s + 9/3 2 / s )]; 
cs - -is [c./.0(4k + 8/3k + 4/3 2 «: + £ - 10/3£ + /3 2 

-c„(-(3 + 4/3 + 9/3 2 )(l + /.)£ + 4(2/3 + 6/3 2 + 3/ s + 2/3/, + 3/3 2 / s )«)]; 
c 9 = -^g[3c s - c„(-3 + 3/3 2 - /, + 2/3/ s )]. 

For S° at WOi: 

ci = -i(-l + /3)((l+/3)c„-2c s ); 

c 2 = -^(-1 + /3)[(-l + /3)c„ + 18(1 + /3)c s / s # 

c 3 = ^(-1 + /3)[c„((-6 - 6/3)/ s -1 + 0)- 6c s ((2 + 2/3)/ s - 1 + /3)]; 

c 4 = £(-1 + /3)[(-l + /3)c u (14 K - 130 - 18(1 + P)c s f s <j>{2n + 0]; 

c 5 = -^[(3 + 2(3 + 3/3 2 ) Cu - (3 + 4/3 + 9/3 2 )c s / s </»(l + /,)]; 

c 6 = i (-1 + /3)[(-l + (3)c u - 18(1 + f3)c s f s ct>}; 

c 7 = 4W(H + 3/3 2 )/ s + 5 - 2/3 + 5/3 2 ) - c s / s ((3 + 4/3 + 9/3 2 )/ s - 3 + 10/3 + 9/3 2 )]; 
c » = -ia4[c„/ s (4K + 8/3« + 4/3 2 k + £ - 10/3^ + /3 2 

-c 8 / 8 #-(3 + 4/3 + 9/3 2 )(l + f s )C + 4((2/3 + 6/3 2 )/ s + 3 + 2/3 + 3/3 2 )«)]; 
c 9 = -^[3c M - c s ((-3 + 3/3 2 )/ s - 1 + 2/3)]. 

For A at WOi: 

ci = -i(-l + /3)((c d + c„)(l + 2/3) - (5 + /3)c s ); 

c 2 = -M-l + P)((7 + H/3)(c d + c u ) + (37 + 35/3)c s / s </>); 

c 3 = -^(-1 + /?)((c + c„)(14 + (8(22 - 3/ s ) + 3/ s ) + c s (-30 - 6/3 + 35/ s + 37/?/,)); 

c 4 = 2M (-1 + /3)((c d + c)(2(-23 + 5/3)k + (17 - 35/3)0 - c s / s </>(58k + 86/3k + 49£ + 23/30); 

c 5 = £((-1 + /3 2 )(c d + c„)(l + /.) + 2(5 + 6/3 + 13/3 2 )c s / s </»); 

c 6 = -^((-1 + /3)((11 + 7/3)(c d + Cu) + (35 + 37/3)c s / s 0)); 

c 8 = ^(2^/^(4(7 + 10/3 + 19/3 2 ) K - (11 + 26/3 + 35/3 2 )0 

+ (-1 + /3)(c d + c)(4(-4 - 2/3 + 5/ s + 7/3/> - (-5 + 11/3)(1 + / s )0); 

c f = skit' 1 + ^)( 2 ( 5 + 0) c « + ( c <* + c «)( 5 + P + n /» + IWs)))- 
At the structure WO2, the sum rules can be expressed in the following form 

Cim SJ L- 8 / 9 ^i(w;)M 2 + c 2 xai" 16/27 ^i(w)M 2 + (c 3 +C4)aEb(iu) + c 5 m sX a 2 L- 16 / 27 E (w) + c 6 m 2 aL- 4 / 9 
+c 7 xabL~ 16/27 + (c 8 + c 9 )m s a 2 ■^+c w ab^ s: 

For £+ at W0 2 : 

c 1 = -i(-l + /3)[(l + /3)c 5 -2c„]; 
c 2 = -|(-l + /3) 2 c s / s </); 

c 3 = -|(-1 + /3)[c s (2 + 2/3 - f a + 0f.) + 2c u (-2 - 2/3 - / s + /3/ s )]; 

c 4 = ^(-1 + /3)[c s / s </> + 2(1 + /3)c„](2« - 0; 

c 5 = -|(-1 + /3)[2c s / s <^ + c(-l -0-f. + /3/ s )]; 

c 6 - £(-1 + /3)[7(1 + /3)c s + c u (-l - /3 - 3/ s + 3/3/ s )]; 

c 7 = - T l 3 (-l + /3) 2 Cs /^; 

cs - 4[c a (-3 + 3/3 2 + 5/ s - 2/3/ s + 5/3 2 / s ) - c„(3 + 4/3 + 9/3 2 - 3/ s + 10/3/ s + 9/3 2 / s )]; 
c 9 = ^[c s / s 0((l - 10/3 + /3 2 )£ + 4(1 + P) 2 k) 

-c„(4(2/3 + 6/3 2 + 3/ s + 2/3/ s + 3/3 2 / s ) K - (3 + 4/3 + 9/3 2 )(l + / s )0]; 

cio = -5gs(-l + WO- + Z 3 )^ - c «(- 3 - 3 ^ - /- + 
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For S° at W0 2 : 

ci = -±(-l + 0)[(l+0)c u -2c s ]; 
c 2 = -±(-l + /3) 2 c M ; 

c 3 = + /9)[cu/.((2 + 20) f a -1+0) + 2c s / s (-(2 + 2)8)/. - 1 + /3)]; 

c 4 = i(-l + /3)[c s / s + 2(1 + /3)c u ](2k - 0; 

c 5 = - |(-1 + /3)[2c u / s + c 5 / s ^»(-(l + /3)/ s - 1 + /3)]; 

c 6 = ^(-1 + /3)[7(1 + /3)c„/ s + c s (-(l + /3)/ s - 3 + 3/3)]; (47) 

c 8 = s[c u / s ((-3 + 3/3 2 )/ s + 5 - 2/3 + 5/3 2 ) - c s / s ((3 + 4/3 + 9/3 2 )/ s - 3 + 10/3 + 9/3 2 )]; 
eg = i[c u / s ((l - 10/3 + /3 2 )£ + 4(1 + /3) 2 K ) 

-c s / s (4((2/3 + 6/3 2 )/ s + 3 + 2/3 + 3/3 2 )k - (3 + 4/3 + 9/3 2 )(l + /.)$]; 
cio = -akC-l + + /5)c„ - c s (-(3 + 3/3)/ s - 1 + 0)}. 

For A at W0 2 : 

ci = + /3)((c d + c u )(l + 2/3) - (5 + 0)c s ); 

C2 = -mi' 1 + /3) 2 ( 2 ( c <i + c u) - c s / s 0); 

c 3 = hC- 1 + /3)((Cd + c «)(2 + /8(-6 + fa) - fa) + +c s (10 + 3/ s + 0(2 + 5/ s ))); 
c 4 = ^(-1 + /?)((- 1 + 3/3) (cj + cu) + (5 + 3/3)c s / s </»)(2 K - £); 
c 5 = M-l + 0)((c d + c M )(-l + 0(-l + /.) - /.) + 4c s / s 0); 

c 6 = g(-l + /3)(-2(5 + 0)c s + {c d + c„)(3 + 7/3 + 5/ s + 3/3/,)); (48) 
c 7 = -skC- 1 + /3) 2 (2(crf + c u ) - c s f s <f)- 

c « = Mi- 1 + P)( c d + c «)(- ? + 0-fs + 13/3/.) - 2c.(-15 + 17/, + 2/3(6 + 7/ s ) + /3 2 (3 + 41/.))); 
c 9 = g| 3 (2c s / s 0(-4(7 + 10/3 + 19/3 2 )k + (11 + 26/3 + 35/3 2 )£) 

-(-1 + 0)(c d + c„)(4(-4 - 2/3 + 5/ s + 7/J/,)« - (-5 + 11/3)(1 + /,)£)); 
c io - -shii' 1 + l 3 )^ + P) c s + (cd + c u )(b + + 11/, + 13/3/,)). 



It is not necessary to list all the coefficients for all 8 
members since the coefficients for other members of the 
octet family can be obtained by appropriate replacements 
of quark contents in the following way: 

• for proton p, replace s quark by d quark in S + . 

• for neutron n, exchange d quark with u quark in 
proton p, 

• for S _ , replace u quark by d quark in E + , 

• for replace u quark by d quark in 5°, 

• for S , change corresponding factors in H . 

Here the conversions between u and d quarks are achieved 
by simply switching their charge factors e u and e<j- More- 
over, The conversions from s quark to u or d quarks in- 
volve setting m s = 0,/ = </i = l,in addition to the 
switching of charge factors. This way, they can also pro- 
vide additional checks for our calculations. 

At this point, we can make some comparisons with pre- 
vious calculations in Ref. 0,0,0]. First, we use general 
interpolating fields where we can vary to achieve the 
best match in sum rules. The previous calculations cor- 
respond to a fixed value of 0=-l. This effect was studied 
in detail in Ref. [l9[, and it was found that = —1.2 is 
the optimal value. Most of our results are at = —1.2. 
Second, we have checked that our sum rules agree with 



I 

those in the previous calculations for the most part. For 
example, for the proton at WEi , we completely agree ex- 
cept for the k— 2£ term in Eq. (2.16) in Ref. y|. In all of 
our sum rules, we have the combination 2k — £ instead of 
k — 2£. For the strange members (£, S, A), they have 8 
terms in the OPE, while we have 10 terms. For the sum 
rules at structure WOi, they have only 3 terms, while we 
have 9 terms. For the sum rules at structure WO2, they 
have only 4 terms, while we have 10 terms. Third, they 
only analyzed the sum rules at structure WEi, while we 
will examine all the structures. Fourth, we use a com- 
pletely different analysis method. 

Before going into the analysis, we would like to point 
out some relations among the correlation functions (or 
OPE) based on symmetries, which lead to the same re- 
lations in the magnetic moments. In exact SU(3)-flavor 
symmetry, it is known that the magnetic moments of the 
octet family are related by (see, for example, Ref. |20j]) 

Ms+ = Mp> 
2/iA = fJ-n, 

A*s- + = -Hp, (49) 
Me- = Ms-, 

fl E o = fl n . 

These relations are borne out in the OPE of our sum rules 
if SU(3)-flavor symmetry is enforced. They are only ap- 
proximately true since SU(3)-flavor symmetry is broken 
by the strange quark. Here we have the advantage of 
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studying the symmetry-breaking effects since the terms 
are explicit in our QCD sum rules. 



IV. SUM RULE ANALYSIS 

The sum rules for magnetic moments have the generic 
form of OPE - ESC = Pole + Transition, or 



U mag (QCD, 0, w, M 2 ) = X 2 N + A 



-Ml/M 2 



(50) 

where QCD represents all the QCD input parameters. 
The task then becomes: given the function IT mas with 
known QCD input parameters and the ability to vary 
/3, find the phenomcnological parameters (magnetic mo- 
ment /xjv, transition strength A, coupling strength Afy, 
and continuum threshold w) by matching the two sides 
over some region in the Borel mass M. A x 2 minimiza- 
tion is best suited for this purpose. It turns out that 
there are too many fit parameters for this procedure to 
be successful in general. To alleviate the situation, we 
employ the corresponding mass sum rules which have a 
similar generic form of OPE - ESC = Pole, or 



ILmass(QCD,l3, W U M 2 



X 2 N e- M «' M 



(51) 



which shares some of the common parameters. Note that 
the continuum threshold may not be the same in the two 
sum rules. By taking the ratio of the two equations, we 
are left with 



n mag (QCD,P, Wl M 2 ) 
U mass (QCD,(3, Wl ,M^ M 2 



This is the form we are going to implement. By plot- 
ting the two sides as a function of 1/M 2 , the slope will 
be the magnetic moment and the intercept the transition 
strength. The linearity (or deviation from it) of the left- 
hand side gives an indication of OPE convergence and 
the role of excited states. The two sides are expected to 
match for a good sum rule. This way of matching the 
sum rules has two advantages. First, the slope, which 
is the magnetic moment of interest, is usually better de- 
termined than the intercept. Second, by allowing the 
possibility of different continuum thresholds, we ensure 
that both sum rules stay in their valid regimes. 

We use the chiral-even mass sum rules in Ref. [2~i| 
which are listed here in the same notation, 

piL- 4 / 9 £; 3 (wi)M 6 +p 2 bL- 4 / 9 E 1 {w 1 )M 2 + p 3 m s aL 4 / 9 
+p 4 a 2 L i / 9 + p b a 2 k v L^ + p 6 m 2 a 2 L~ 2 / 27 ^ 



\ 2 N e- M »/ M 



The coefficients for N are: 

Pi-m(5 + /? + 5/3 2 ); 
P2= 2ie(5 + /3 + 5/3 2 ); 

P4 = £(7-2/3-5/3); 
P5 = 0; 

p 6 --^(13-2/3-ll/3 2 ). 



M- 



(53) 



(54) 



For A: 

p 1 = ^(5 + 2/? + 5/? 2 ); 
P2= 2le (5 + 2/3 + 5/3 2 ); 

P3 = £((20 - 15/.) - (16 + 6/.)/3 - (4 + 15/ s )/3 2 ); 
P4 = i((4/ s - 5 - 6i) + (4 + if s )(3 + (4/, + 1 + 6t)/3 2 ); 
P5 = MlOfs + 11) + (2 - 8/.)/3 ~ (2/ s + 13)/? 2 ); 
Pe = m((- l6 fs - 23) + (8/ s - 2)(3 + (8f s + 25)/3 2 ). 

(55) 

For E: 



Pi 

Pi 



M(5 



2/3 + 5/J 2 ); 
f 2/3 + 5/3 2 ); 



64 * 

256 ( 5 ' 

p 3 = g((12 - 5/.) - 2f s (3 - (12 + 5/ s )/? 2 ); 
P4 = -h(( 4 /- + 21 + 18*) + 4/ s /3 + (4/, - 21 - 18t)P 2 ); 
P5= 2l((6/ s + l)-2/5-(6/ s -l)/3 2 ); 
p 6 = -i((12/ 5 + 1) - 2/3 - (12/, - l)/3 2 ). 

(56) 

For S: 



Pi 

P2 
P3 
Pi 
P-> 
P6 



6|( 5 ^ 
356 ( 5 ' 

i((2- 

-m((15' 



2/? + 5/? 2 ); 
-2/3 + 5/3 2 ); 

-/,)-2/ s /3-(2 + / s )/3 2 ); 

f s + 18t) - 10/./3 - (15 + f s + 18t)/3 2 ); 
= ^/ s ((/ s + 6)-2/ s/ 3+(/ s -6)/3 2 ); 
= -&f,((f. + 12) ~ 2fsP + (fs ~ 12)P 2 ). 

(57) 

M 2 

The function t is defined as t = In — Jem with Jem « 
0.577 the Euler-Mascheroni constant. Use this mass sum 
rule for the coupling constant for magnetic moment sum 
rules, and move to the left side, then only the magnetic 
moment pn and the transition strength A dependence 
are on the right hand side. 

We use the Monte-Carlo procedure first introduced in 
Ref. [HI to carry out the search which allows a rigorous 
error analysis. In this method, the entire phase-space of 
the input QCD parameters is explored simultaneously, 
and is mapped into uncertainties in the phenomcnolog- 
ical parameters. This lead to more realistic uncertainty 
estimates than traditional approaches. 

First, a set of randomly-selected, Gaussianly- 
distributed condensates are generated with a assigned 
uncertainties. Here we give 10% for the uncertainties 
of input parameters, and this number can be adjusted 
to test the sensitivity of the QCD parameters. Then 
the OPE is constructed in the Borel window with evenly 
distributed points Mj. Note that the uncertainties in 
the OPE are not uniform throughout the Borel window. 
They are larger at the lower end where uncertainties in 
the higher-dimensional condensates dominate. Thus, it 
is crucial that the appropriate weight is used in the cal- 
culation of \ 2 . For the OPE obtained from the k'th set 
of QCD parameters, the x 2 P er degree of freedom is 

X\ ^ [Tir E (M 2 ,(3,w, Wl ) -n^(M 2 ,^A)] 2 



iV, 



DF 



E 



(n B - n p ) a 2 PE (Mj) 



(58) 

where Ii OPE refers to the LHS of Eq. {52]) and W hm its 
RHS. The integer n p is the number of phenomenologi- 
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cal search parameters. In this work, tib=51 points were 
used along the Borel axis. The procedure is repeated for 
many QCD parameter sets, resulting in distributions for 
phenomenological fit parameters, from which errors are 
derived. In practice, 200 configurations are sufficient for 
getting stable uncertainties. We used about 2000 sets 
to resolve more subtle correlations among the QCD pa- 
rameters and the phenomenological fit parameters. This 
means that each sum rule is fitted 2000 times to arrive 
at the final results. 

The QCD input parameters are given as follows. 
The condensates are taken as a = 0.52 GeV 3 ,b — 
1.2 GeV 4 ,ml = 0.72 GeV 2 . For the factorization vio- 
lation parameter, we use k v — 2.0. The QCD scale pa- 
rameter is restricted to Aqcd — 0.15 GeV. The vacuum 
susceptibilities have been estimated in studies of nucleon 
magnetic moments [3J, [5j, |9|, [12J , but the values vary de- 
pending on the method used. We use \ = —6.0 GeV' 2 
and k = 0.75, £ = —1.5. Note that x is almost an order of 
magnitude larger than k and £, and is the most important 
of the three. The strange quark parameters are placed 
at m s = 0.15 GeV, f = 0.83, <p = 0.60 d ©]. These 
input parameters are just central values. We will explore 
sensitivity to these parameters by assigning uncertainties 
to them in the Monte-Carlo analysis. 



V. RESULT AND DISCUSSION 

We have 24 sum rules in total to analyze: 3 for each 
member of the octet. For each sum rule, we have in prin- 
ciple 5 parameters to determine: fi, A, w, w\, f3. But a 
search treating all five parameters as free does not work 
because there is not enough information in the OPE. In 
fact, the freedom to vary (3 can be used as an advantage 
to yield the optimal match. We find that (3 = —1.2 gives 
the best match in most cases. This agrees with the value 
suggested in Ref. [19(. One exception is the proton: we 
found a better solution at (3 = —0.4 than at f3 = —1.2. 
Another parameter that can be used to our advantage is 
the continuum threshold w\ for the corresponding mass 
sum rule. We fix it to the value that gives the best solu- 
tion to the mass sum rule independently. The following 
values for w\ are used: for the nucleon, w± = 1.44 GeV; 
for A, wi = 1.60 GeV; for E, w t = 1.66 GeV; for S, 
wi = 1.82 GeV. In this way the magnetic moment sum 
rule and the mass sum rule can stay in their respective 
valid Borel regimes. This leaves us with three parame- 
ters: /i, A, w. Unfortunately, a three-parameter search is 
either unstable or returns values for w smaller than the 
particle mass, an unphysical situation. Again we think 
this is a symptom of insufficient information in the OPE. 
So we are forced to fix the continuum threshold w that 
corresponds to the best match for the central values of 
the QCD parameters. 




2.5 3.0 




Transition Amplitude A 

FIG. 4: Histogram for the proton magnetic moment (top) and 
transition amplitude (bottom) obtained from Monte-Carlo fits 
of Eq. [32j) at WEj for 2000 QCD parameter sets. They are 
based on 10% uncertainty given to all the QCD input param- 
eters. 



A. The Sum Rule at WEi 

The results determined this way at the WEi structure 
are displayed in Table HI The Borel window is deter- 
mined by the following two criteria: OPE convergence 
which gives the lower bound, and ground-state domi- 
nance which gives the upper bound. It is done iteratively. 
For each value of /3, we adjust the Borel window until the 
best solution is found. We see that our calculated mag- 
netic moments agree with experiment fairly well within 
error bars. 

We stress that the errors are derived from Monte-Carlo 
distributions which give the most realistic estimation of 
the uncertainties. An example of such distributions is 
given in Fig. 2] We see that they are roughly Gaussian 
distributions. The central value is taken as the average, 
and the error is one standard deviation of the distribu- 
tion. We found about 10% accuracy for the magnetic mo- 
ments in our Monte-Carlo analysis, resulting from 10% 
uniform uncertainty in all the QCD input parameters. 
Of course, the uncertainties in the QCD parameters can 
be non-uniform. For example, we tried the uncertaint 



assignments (which are quite conservative) in Ref. 
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TABLE I: Results for the magnetic moment of octet baryons from the QCD sum rule in Eq. (|32[) (structure WEi). The seven 
columns correspond to, from left to right: particle, (3 value, Borel region in which the two sides of the QCD sum rule are matched, 
continuum threshold, transition strength, extracted magnetic moment in unit of nuclear magnetons, and experimental value. 
The errors are derived from 2000 samples in the Monte-Carlo analysis with 10% uncertainty on all QCD input parameters. 





ft 
P 


xxe^ioii 


w 


A 


fl B 


Jjjxp. 






l bev ) 


Oev ) 


((s.e\ ) 


In \ 




V 


-1.2 


0.7 to 0.9 


1.40 


1.46 ± 0.34 


3.01 ± 0.24 


2.79 




-0.4 


0.8 to 1.2 


1.60 


0.74 ± 0.13 


2.82 ± 0.26 


2.79 


11 


-1.2 


0.7 to 1.1 


1.40 


-0.19 ± 0.09 


-1.97 ± 0.15 


-1.91 


A 


-1.2 


1.1 to 1.2 


1.60 


-0.45 ± 0.05 


-0.56 ± 0.15 


-0.61 


E+ 


-1.2 


1.1 to 1.3 


1.85 


0.56 ± 0.06 


2.31 ± 0.25 


2.45 


E° 


-1.2 


1.0 to 1.6 


1.80 


0.07 ± 0.02 


0.69 ± 0.07 


0.65 


£" 


-1.2 


1.3 to 1.8 


1.80 


-0.08 ± 0.01 


-1.16 ± 0.10 


-1.16 




-1.2 


1.6 to 1.9 


2.15 


-0.16 ± 0.01 


-1.15 ± 0.05 


-1.25 




-1.2 


1.1 to 1.4 


2.00 


-0.34 ± 0.02 


-0.64 ± 0.06 


-0.65 



and found about 30% uncertainties in our output. 

To gain a better appreciation on how the QCD sum 
rules produce the results, we show Fig. O using the pro- 
ton as an example. There are three graphs in this figure 
to give three different aspects of the analysis. The first 
graph shows how the two sides of Eq. (|52|) match over 
the Borel window, which should be linear as a function 
of 1/M 2 according to the right-hand side of this equa- 
tion. Indeed, we observe excellent linear behavior from 
the OPE side (LHS). The match is almost perfect (barely 
distinguishable between the solid and dotted lines). The 
slope gives the magnetic moment /i, and the intercept 
give the transition contribution A. We find that the in- 
clusion of A is important in producing the best match. 
Also plotted arc the individual contributions from u and 
d quarks. We see that for the proton, the u-quark con- 
tribution is the dominant one, which is expected because 
it is doubled-represented in the proton (uud). We define 
the slope from an individual quark contribution as the 
effective magnetic moment of that quark in the particle. 

The second graph in Fig. [3] shows how the various 
terms in the OPE contribute to determination of mag- 
netic moments. The M° term, which contains the con- 
tributions from the condensates x® 2 an d b, plays an im- 
portant role. It is the leading contribution in the region 
below M 2 < 1.2 GeV 2 . For this reason, the sum rule at 
WEi is expected to have good spectral properties. In- 
deed this is confirmed in the third graph where we plot 
the three terms in the phenomenological side (pole, tran- 
sition, and excited) as a function of M 2 . The ground- 
state pole is dominant (over 70% of the RHS at the low 
end of the Borel window). The excited-state contribution 
starts small, then grows with M 2 , as expected from the 
continuum model. The transition contribution is small in 
this sum rule. It is consistently smaller than the excited- 
state contribution and has a weak dependence on the 
Borel mass. 



B. The Sum Rule at WOi 

Next, we analyze the sum rule in Eq. (|41 [) at the struc- 
ture WOi, using the same procedure. Table HT1 displays 
the results extracted from this sum rule. The magnetic 
moments have larger errors than those from WEi: about 
15% as opposed to 10%. The agreement with experiment 
is reasonable (with the exception of £~), but not as good 
as those from WEi . We had to search a wider region in (3 
to find the best match. The transition contribution (A) 
are larger for the strange particles, as well as their errors. 

Fig. [5] shows the details of the analysis in the case of 
the proton. The matching is very good, as indicated in 
the top graph. The middle graphs shows that the lead- 
ing contribution in the OPE (M 4 term) is x a i followed 
by the quark condensate a (M 2 term). The condensate 
Xab (M° term) and ab (M~ 2 term) are very small in this 
sum rule. The bottom graph reveals a surprising result: 
the excited-state dominates over the pole and the tran- 
sition. As a result, this sum rule is less reliable. This 
is the reason why the results from this sum rule are not 
as good as those from WEi. This sum rule also shows 
the importance of checking the individual terms in the 
phenomenological side, in addition to looking at the best 
match of the two sides. In this case, there is no pole domi- 
nance, even though the leading terms is non-perturbative 
and the match is almost perfect. 

C. The Sum Rule at WO z 

Finally, we present the results from the sum rule in 
Eq. gS]) at the structure W0 2 in Table HO The agree- 
ment with experiment is not as good as the other two 
sum rules. For example, E~ and S~ have the wrong 
sign. Fig. [7] shows the details of the analysis for the pro- 
ton. The matching is very good, as indicated in the top 
graph. The middle graph shows that the leading contri- 
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TABLE II: Similar to Table|U but for the QCD sum rule in Eq. |gT]| (structure WOi). 





P 


Region 


w 


A 




Exp. 






(GeV) 


(GeV) 


(GeV" 1 ) 




(/x N ) 


p 


-0.8 


1.4 to 1.6 


1.50 


1.03±0.16 


2.67 ±0.16 


2.79 


n 


-1.2 


1.2 to 1.4 


1.40 


-0.54 ±0.11 


-1.70 ±0.11 


-1.91 


A 


-1.2 


1.0 to 1.3 


1.60 


-2.3 ±0.31 


-0.62 ±0.17 


-0.61 


E+ 


-0.6 


1.5 to 1.7 


1.60 


1.60 ±0.50 


2.46 ±0.40 


2.45 


E° 


-1.0 


1.4 to 1.6 


1.70 


0.23 ±0.11 


0.58 ±0.09 


0.65 


S" 


-0.6 


1.2 to 1.5 


1.60 


-2.04 ±0.34 


-0.57 ±0.26 


-1.16 




-1.2 


1.2 to 1.4 


2.10 


-1.52 ±0.21 


-1.27 ±0.15 


-1.25 




-0.2 


1.7 to 1.8 


1.90 


-0.4 ±0.41 


-0.49 ±0.29 


-0.65 



bution in the OPE is M 2 with a coefficient of xa, followed 
by the M° term. The 1/M 2 term is slightly negative, 
while the 1/M 4 term is very small. The bottom graph 
shows that the excited-state dominates over the pole and 
the transition, like the sum rule from WOi, but the rel- 
ative size of the pole is much larger. Since the WO2 sum 
rule has power corrections up to 1/M 4 , it is expected to 
be more reliable than the WOi sum rule. But our anal- 
ysis shows that this advantage is offset by the smallness 
of the 1/M 2 and 1/M 4 terms. As a result, the reliability 
of the WO2 sum rule is about the same as the WOi sum 
rule. 

We have performed the same analysis for all the mem- 
bers and all three structures. Fig. [5] shows the graphs for 
the neutron. In this case, the slope is negative. Again, 
the sum rule at WEi has excellent convergence proper- 
ties. The WOi sum rule has a good match, but the pole 
is less than the excited-state. The WO2 sum rule does 
not have a good match. Fig.[9]shows the case for the E , 
which has all three quark contributions (u, d and s). 

Overall, based on the quality of the match, the broad- 
ness of the Borel window and its reach into the lower end, 
the size of the continuum contribution, and the OPE con- 
vergence, we find that the sum rule at WEi is the most 
reliable of the three sum rules. 



D. Some physics discussions 

Based on the results of our comprehensive analysis, we 
conclude that the QCD sum rules at WEi in Eq. ([5]) 
are the most reliable. Here we discuss some physics im- 
plications of the results extracted from them (Table [T| . 
First we look at some ratios of magnetic moments listed 
in Table HVl Here we compare our results with those 
from the SU(6) symmetry, lattice calculations [22j and 
experiment. From the table, we see that the QCD sum 
rule results compare well against other approaches and 
experiment. They agree a little better with experiment 
than the lattice results. Furthermore, our QCD rum rule 
results are an improvement over the ratios from previous 
QCD sum rule calculations in Ref. []|. 



Next, we consider a few sum rules among magnetic 
moments that have been discussed in literature. They 
reveal interesting quark dynamics in the baryons. They 
are mostly based on SU(6)-symmetry considerations in 
the quark model. We begin with the sum rule [23| 



3(p - ^ 



(p + 3A) 



P 



(59) 



It assumed 'baryon independence' of quark moments: the 
independence of which baryon the same quark is in, a 
concept first mentioned by Franklin 24j. In other words, 
each quark is not sensitive to the environment it resides in 
and quarks in different spin states had the same effective 
moments. In this sum rule, it also assumed that the s 
quark moment in the A is the same as in the E and S, 
even though the spin states are different. This rum rule 
is violated the most (by a factor of 5) using present values 
of magnetic moments. This large violation is mostly due 
to the small difference in the denominator S~ — S° which 
magnifies the apparent discrepancy. So this is not a good 
way of testing baryon independence. Our determination 
of the ratio of the left hand side over the right hand 
side (LHS/RHS) is 7.4(1.0), compared to the experiment 
value of 5.7(4) and the lattice calculation of 4.1(1.5). The 
error are added in quadrature in forming these ratios. 

One interesting sum rule [29| which is fairly accurately 
satisfied is 



p + n = 3A± -(E H 



(SP+S-). 



(60) 



It is derived using SU(3)-flavor symmetry to characterize 
the A wave function Our determination of the LHS /RHS 
ratio is 1.24(79), compared to the experimental values of 
1.22(5) and the lattice calculation of 1.45(27). 
Another sum rule, first derived by Franklin [3], 



p — n = E — E 



(61) 



is another test of baryon independence of the quark mo- 
ments. The strange quarks approximately cancel, leav- 
ing only u and d quarks. We give the ratio of 1.20(0.68), 
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TABLE III: Similar to Table |U but for the QCD 


sum rule in Eq. 


(|45l) (structure WO 


2). 




P 


Region 


w 


A 


Hb 


Exp. 






(GeV) 


(GeV) 


(GeV -2 ) 


(m«) 


(M«) 


p 


-1.2 


1.3 to 1.4 


1.60 


0.42 ±0.07 


2.53 ±0.13 


2.79 


n 


-1.2 


1.7 to 1.9 


1.60 


0.26 ±0.03 


-1.73 ±0.28 


-1.91 


A 


-1.2 


1.3 to 1.4 


1.90 


-0.85 ±0.22 


-0.53 ±0.17 


-0.61 


E+ 


-1.0 


1.2 to 1.4 


1.70 


0.21 ±0.05 


1.50 ±0.11 


2.45 


E° 


-1.2 


1.0 to 1.6 


1.60 


-0.02 ±0.02 


1.04 ±0.07 


0.65 


£- 


-1.2 


1.4 to 1.6 


1.65 


0.0 ±0.01 


0.61 ±0.06 


-1.16 




-1.2 


1.6 to 2.0 


2.1 


-0.08 ±0.03 


-1.20 ±0.23 


-1.25 




-1.2 


1.4 to 1.5 


2.2 


0.22 ±0.09 


0.80 ±0.21 


-0.65 



while the experiment measurements gives 1.133(10) and 
lattice results is 1.07(8) for this ratio. 
Finally, the Sachs sum rule [26[ 

3(p + n) = - £~ - S~ ±S° (62) 

is satisfied by the more general extension of SU(6) sym- 
metry and is another test of baryon independence of 
quark moments. This sum rule is just the sum of the 
two separate sum rules proposed earlier by Franklin [25| . 
Our result gives the ratio of 0.86(35) which is in agree- 
ment with the experiment measurements of 0.881(11). 
Whereas it yields an opposite violation with a ratio of 
1.29(20) from lattice moments. This contradiction is pos- 
sibly because the Sachs sum rule may be sensitive to dy- 
namics not included in the lattice calculation. 



E. Individual quark contributions 

To gain a deeper understanding of the dynamics, it is 
useful to consider the individual quark sector contribu- 
tions to the magnetic moment. In our approach, we can 
easily dial individual quark contributions to the QCD 
sum rules. For example, to turn off all u-quark contri- 
butions, we set the charge factor c u = 0. To turn off all 
s-quark contributions, we set c s — 0, m s = 0, / = 1, 
and <f> = 1. We can extract a number corresponding to 
each quark contribution from the slope of Eq. (|S"2"|) as a 



TABLE IV: A comparison of selected ratios of magnetic mo- 
ments. 



Ratio 


Sum Rule 


SU(6) 


Lattice 


Exp. 




result 


symmetry 


results 




n/p 


-0.70(9) 


-2/3 


-0.63(5) 


-0.68 


ST/E+ 


-0.50(7) 


-1/3 


-0.37(3) 


-0.47 




0.56(6) 


1/2 


0.58(5) 


0.52 


H°/A 


2.05(56) 


2 


2.4(5) 


2.04 


B-/A 


1.14(32) 


1 


1.0(2) 


1.13 



function of \/M 2 . We call this the raw individual quark 
contributions to the magnetic moments. 



Table [V] gives the result of raw individual u, d and 
s quark sector contributions to the magnetic moments 
from the QCD sum rules at WEi. It is compared with 
the lattice QCD result in [2^]. The lattice results were 
rescaled later in (2?| and we use the rescaled results for 
the comparison. Our results agree with the lattice results 
reasonably well. The biggest discrepancy is that our light 
quark moments in S° and S~ are small. In £ , and A 
the total of u and d quark moments agree very well with 
lattice data for light quark moments. Note that only a 
combined number for light quark (fxi = fx u + /J,d) was 
given on the lattice [22J . To compare with our separated 
u and d quark numbers, we break up the lattice number 
by using the relation fj, u = — 2/x^. 



In the simple quark model [281 ], the magnetic moment 
of the proton is given by p — |/i u — ^fid- In the SU(2) 
limit, there is the relation fi u = —2^. Our results for 
S° and A, the two uds structure baryons, agree with 
this very well. This relation also suggests that the ra- 
tio of the u quark and d quark contributions in pro- 
ton is = 8. Our QCD sum rule ratio is 
10.6 for (3 = —1.2, which agrees with the lattice ratio 
of 10.3(7) at K\. For neutron the quark model ratio is 
^d/{— \nu) = 2, and our result is 3.64(40) for this ra- 
tio, which differs from the SU(6) spin-flavor symmetry 
prediction. The lattice ratio is 2.6. The bigger ratio rep- 
resents an enhancement of the doubly represented quark 
contribution. These results suggest that although the to- 
tal magnetic moments ratios agree with SU(6) ratios (see 
Table lTV) , the underlying quark dynamics are really quite 
different from different individual quark contributions. 



To include the possible nonstatic effect, Franklin [2^| 
proposed a generalization of SU(6) results so that non- 
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TABLE V: Individual quark contributions and total magnetic moments in unit of nuclear magnetons extracted from the QCD 
sum rules at WEi (denoted by the subscript SR) are compared with those from a lattice QCD calculation (denoted by the 
subscript LAT) 





PSR 


Mlat 


,, d 

PSR 


,, d 
Plat 


,,s 
PSR 


Mlat 


Total 


, Total 
PLAT 


p 


2.46(25) 


2.59(24) 


0.38(4) 


0.15(9) 








2.82(26) 


2.79 


It 


-0.42(9) 


-0.31(20) 


-1.53(13) 


-1.32(14) 








-1.97(15) 


-1.60(21) 


A 


0.18(7) 


0.09(6) 


-0.09(4) 


-0.03(2) 


-0.63(15) 


-0.54(7) 


-0.56(15) 


-0.50(7) 


£+ 


2.06(25) 


2.33(29) 








0.15(2) 


0.08(5) 


2.31(25) 


2.37(18) 


E° 


1.09(12) 


0.87(9) 


-0.54(6) 


-0.29(3) 


0.14(2) 


0.08(5) 


0.69(7) 


0.65(6) 


£- 








-1.25(10) 


-1.14(14) 


0.10(2) 


0.08(5) 


-1.16(10) 


-1.08(10) 




0.02(2) 


-0.44(6) 








-1.08(7) 


-0.73(6) 


-1.15(5) 


-1.17(9) 










0.03(2) 


0.22(3) 


-0.67(6) 


-0.73(6) 


-0.64(6) 


-0.51(7) 



static components are the same for each octet baryon: 



(63) 



4 1 / 

P — 3 fJ~u 3 Md' 
n = - §/4, 
y+ - ±n - in' 

v- — 4 1 ,,' 

^ — 3 Pd 3^8! 

"° - — II - ±u' 

a = 3^ s — 3^d' 

A = 



where quark symbols refer to quark moment contribu- 
tions including nonstatic effect. The magnetic moment 
contribution of the unlike quark in the baryon is primed. 

In order to compare with the effective quark moments 
defined in the quark model, we convert our raw quark 
moments in Table quark in a similar fashion. Take the 
proton for example, we define effective quark moments 
/i u and \i! d from QCD sum rules by usr = ^fi u , dsn = 
\fi' d - One can define effective moments for quarks in 
other baryons in a similar manner. 

In the proton, fi u is 1.84(18) while in £+ 1.54(18) and 
in E° is 1.63(18). The magnetic moments list here are 
all in unit of (in. In the neutron, /i' u is 1.25(27), which 
is smaller than fi u . This is an example that effective 
quark magnetic moment is sensitive to the environment 
the quark resides in. We find the following relation from 
our results 



(64) 



For n d , in n it is -1.15(11) while for S° is -0.82(8) and 
for £- is -0.93(8). From p we can get \x' d = -1.15(12). 
It has a similar relation for the absolute value 



(65) 



For fi s , it is -0.81(6) in S° and -0.50(5) in E - . For 
fi' s , it is -0.45(5), -0.42(5) and -0.30(5) in £+, E° and 
E~, respectively From A we can get fi!' s about —0.63(15). 

For these individual quark effective moments, we no- 
tice that in general we have the following relation 



fi s < fl d < fl u . 



(66) 



It is expected because of the quark mass effects, which is 
analogous to those seen in the electric properties. 

Another way of looking at the individual effective 
quark moments is by expressing them in terms of baryon 
magnetic moments using Eq. (|63|) and isospin symmetry. 
For example, the d-quark effective moment can be ex- 
pressed as 

Md = -~(2p + n) = -(£--£+). (67) 

Our result indicates fid = —0.92(14) < —0.86(7) using 
the baryon moments in Table [V] or in Table HI It agrees 
well with the experimental moments fid = —0.918 < 
-0.894(7) and the lattice result of \± d = -1.00(5) < 
-0.86(6). 

Similarly the dl quark effective moment is 



Vd=P 



2n = Z°- 5" 



(68) 



We have fi' d = -1.12(39) < -0.51(8). It also agree very 
well with the experimental moments fi' d = —1.003 < 
—0.578(26) although the two sides of this sum don't agree 
well. 

The strange quark can be isolated as 



KH°+2H- 



fi' s = -£+ - 2£- 
tf = A- a. 



(69) 



Our results are compared with experiment and lattice 
calculation in Table IVT1 Not only the results agree with 
the experiment and lattice calculation, they also roughly 
agree with the strange quark moments from individual 
quark moments. 

Now, we can look at the quark moment difference such 
as jjf s — fi' d and fi s — fid- According to simple quark model, 
the difference should be the same and approximately as 
0.36. In fact, from 



^ = 3(p-E H 



(70) 



17 



Pole+Tran 
OPE-ESC 




0.6 O.B 1.0 1.2 1.4 1.8 

l/M 2 




0.6 0.6 1.0 1.2 1.4 1.6 



FIG. 5: Analysis of the QCD sum rule in Eq. ([32]) (structure 
WEi) for the proton at f3 = —0.4 according to Eq. (f52|) . In 
the top, the pole plus transition terms (solid lines) are com- 
pared against the OPE minus the excited-state contributions 
(dashed lines) as a function of l/M 2 (the two should match 
for an ideal sum rule). Also plotted are the individual contri- 
butions from u (long-dashed lines) and d (dot-dashed lines) 
quarks. In the middle, the total in the OPE side and its var- 
ious terms are plotted as a function of M 2 . In the bottom, 
the 3 terms in the phenomenological side: pole (solid), tran- 
sition (long-dashed), and excited (dot-dashed) are plotted as 
a function of M 2 . 



/*.-W=|(S°-n), (71) 

QCD sum rules give 1.53(8) and 0.61(12). The exper- 
iment gives 1.12(7) and 0.495(11), while lattice calcu- 
lation gives 1.3(5) and 0.32(17) respectively. The large 
result of fi' s — fi' d is difficult to reconcile with any sim- 
ple model. And Eq. [70] is a poor way to measure the 




.38 .40 .42 .44 .46 .48 .50 



l/M' 

J— Totil 1 




1.9 2.0 2.1 2.2 2.3 2.4 2.5 2.6 



FIG. 6: Similar to Fig.0 but at structure WOi and j3 = -1.2. 



TABLE VI: The result for s s' and s" quark effective moments 
from QCD sum rule, experiment and lattice calculation. 





Ms 


n 


/ 

Ms 


SR 


-0.61(3) < 


-0.60(15) < 


0.01(32) 


Exp. 


-0.651(12) < 


-0.57(4) < 


-0.107(36) 


LAT. 


-0.55(4) » 


-0.54(4) < 


-0.24(17) 



difference of s and d quark contribution. 



F. Correlations 

Our Monte-Carlo analysis affords the opportunity to 
study the correlations between any two parameters since 
the entire QCD input phase space is mapped into the 
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FIG. 7: Similar to Fig.0 but at structure WQ 2 and (3 = -1.2. 



phenomenological output space. This correlation can be 
explored by a scatter plot of the two parameters of inter- 
est. 

Fig. ITTJI shows the scatter plots for the proton magnetic 
moment at structure WEi . We see that the magnetic mo- 
ment has a strong correlation with the vacuum suscepti- 
bility x- It is a negative correlation meaning larger x ( hi 
absolute terms since x is negative) leads to smaller \xb- 
A slight negative correlation with the mixed condensate 
and a slight positive correlation with another vacuum 
susceptibility k are also observed. Precise determination 
of the QCD parameters, especially for those that have 
strong correlations to the output parameters, is crucial 
for keeping the uncertainties in the spectral parameters 
under control. 

Fig. [11] shows a similar plot for the S° at structure 
WEi . Here we focus on the three vacuum susceptibilities 



and the three parameters that define the strange quark 
(m s , / and </>). The correlations with the other conden- 
sates are similar to the proton and are not shown. A 
negative correlation with x exists, but not as strong as 
that for the proton. A slight positive correlation with m s 
and / is also observed. 



VI. CONCLUSION 

We have carried out a comprehensive study of the mag- 
netic moment of octet baryons using the method of QCD 
sum rules. We derived a new, complete set of QCD sum 
rules using generalized interpolating fields and examined 
them by a Monte-Carlo analysis. Here is a summary of 
our findings. 

We proposed a new way of determining the magnetic 
moments from the slope of straight lines. We find this 
method more robust than from the normalization (inter- 
cept) or from looking for 'flatness' as a function of Borel 
mass. The linearity displayed from the OPE side matches 
almost perfectly with the phenomenological side in most 
cases. The method also demonstrates clearly that the 
transition terms caused by the external field in the in- 
termediate terms cannot be ignored. They are needed to 
make the two sides of a QCD sum rule match. 

Out of the three independent structures, we find that 
the sum rules from the WEi structure are the most re- 
liable based on OPE convergence and ground-state pole 
dominance. The QCD sum rules from this structure are 
in Eq. ([52]); its prediction are found in Table HI and con- 
vergence properties are displayed in Fig. [5] They should 
be considered as the best results in this work. The ex- 
tracted magnetic moments are in good agreement with 
experiment. These results are used to shed light on a 
variety of magnetic moment sum rules based on SU(6) 
spin-flavor symmetries in the quark model, along with 
experiment and lattice QCD. Reasonable results from the 
other two structures (WOi and WO2) are obtained for 
the first time, but they are less reliable due to poor con- 
vergence properties. 

Our Monte-Carlo analysis revealed that there is an un- 
certainty on the level of 10% in the magnetic moments if 
we assign 10% uncertainty in the QCD input parameters. 
It goes up to about 30% if we adopt the conservative 
assignments that have a wide range of uncertainties in 
Ref. [IH . The Monte-Carlo analysis also revealed some 
correlations between the input and output parameters. 
The most sensitive is the vacuum susceptibility x- So 
a better determination of this parameter can help im- 
prove the accuracy on the magnetic moments and other 
quantities computed from the same method. 

We also isolated the individual quark contributions to 
the magnetic moments. These contributions provide in- 
sight into the rich dynamics in the baryons. By compar- 
ing them with the simple quark model and lattice QCD 
results, we reveal the effects of SU(3)-flavor symmetry 
breakings in the strange quark, the environment sensi- 
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dated and improved determination of the magnetic mo- 
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sophistication as the decuplet baryons |9j. One possi- 
ble extension along this line is a calculation of the N 

to A electromagnetic transition amplitudes which to our This work is supported in part by U.S. Department of 
knowledge have not been studied in this method. Energy under grant DE-FG02-95ER-40907. 



[1] M.A. Shifman, A.I. Vainshtein and Z.I. Zakharov, Nucl. 

Phys. B147, 385, 448 (1979). 
[2] I.I. Balitsky and A.V. Yung Phys. Lett. B129, 328 

(1983) . 

[3] B.L. Ioffe and A.V. Smilga, Nucl. Phys. B232, 109 

(1984) . 

[4] B.L. Ioffe and A.V. Smilga, Phys. Lett. B133, 436 
(1983); 

[5] C.B. Chiu, J. Pasupathy, S.L. Wilson, Phys. Rev. D33 , 
1961 (1986). 

[6] J. Pasupathy, J. P. Singh, S.L. Wilson, and C.B. Chiu, 

Phys. Rev. D36 , 1442 (1986). 
[7] S.L. Wilson, J. Pasupathy, C.B. Chiu, Phys. Rev. D36 , 

1451 (1987). 

[8] S. Zhu, W. Hwang and Z. Yang, Phys. Rev. D57, 1527 
(1998). 

[9] F.X. Lee, Phys. Rev. D57, 1801 (1998). 
[10] F.X. Lee, Phys. Lett. B419, 14 (1998). 
[11] J. Dey, M. Dey and A. Iqubal, Phys. Lett. B477, 125 
(2000). 

[12] M. Sinha, A. Iqubal, M. Dey and J. Dey, Phys. Lett. 

B610, 283 (2005). 
[13] A. Samsonov, Phys. Atom. Nucl. 68, 114 (2005). 
[14] T. M. Aliev, I. Kanik and M. Savci Phys. Rev. D68, 



056002 (2003). 

[15] T. M. Aliev, A. Ozpineci and M. Savci, Phys. Rev. D66, 
016002 (2002). 

[16] T.M. Aliev, A. Ozpineci and M. Savci, Phys. Lett. B516, 
299 (2001). 

[17] T.M. Aliev, A. Ozpineci and M. Savci, Nucl. Phys. A678, 
443 (2000). 

[18] B. L. Ioffe, Nucl. Phys. B 188, 317 (1981); Nucl. Phys. 

B 191, 591(E) (1981). 
[19] D. B. Leinweber, Ann. of Phys. (NY.) 254, 328 (1997). 
[20] S. Coleman and S. Glashow, Phys. Rev. Lett. 6, 423 

(1961). 

[21] F. X. Lee and X. Liu, Phys. Rev. D66, 014014 (2002). 

[22] D. B. Leinweber, Phys. Rev. D43, 1659, (1991). 

[23] H. J. Lipkin, Phys. Rev. D24, 1437 (1981). 

[24] J. Franklin, Phys. Rev. 182, 1607 (1969). 

[25] J. Franklin, Phys. Rev. D20, 1742 (1979). 

[26] R. G. Sachs, Phys. Rev. D23, 1148 (1981). 

[27] D. B. Leinweber, Phys. Rev. D45, 252 (1992). 

[28] D. H. Perkins, Introduction to High Energy Physics 

(Cambridge, 4th ed., 2000). 
[29] J. Franklin, Phys. Rev. D29, 2648 (1984). 



20 




21 





-10 -8 -6 -4 -2 0.4 0.6 0.8 1 1.2 -2.5 -2 -1.5 -1 -0.5 

X (GeV 2 ) k | (GeV 2 ) 

FIG. 10: Scatter plots showing correlations between the magnetic moment and the QCD parameters for the proton at structure 
WEi. They are obtained from 2000 Monte-Carlo samples with 10% uncertainty on all the QCD parameters. 




